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Abstract 

The evolution of complex molecular traits such as disulphide bridges often requires 
multiple mutations. The intermediate steps in such evolutionary trajectories are likely 
to be selectively neutral or deleterious. Therefore, large populations and long times 
may be required to evolve such traits. We propose that errors in transcription and 
translation may allow selection for the intermediate mutations if the final trait provides 
a large enough selective advantage. We test this hypothesis using a population based 
model of protein evolution. If an individual acquires one of two mutations needed 
for a novel trait, the second mutation can be introduced into the phenotype due to 
transcription and translation errors. If the novel trait is advantageous enough, the 
allele with only one mutation will spread through the population, even though the gene 
sequence does not yet code for the complete trait. The first mutation then has a higher 
frequency than expected without phenotypic mutations giving the second mutation a 
higher probability of fixation. Thus, errors allow protein sequences to "look-ahead" for 
a more direct path to a complex trait. 
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1 Introduction 



According to a central principle of molecular evolution, the likelihood that a given mutation 
occurs is independent of the mutation's phenotypic consequences. Organisms cannot choose 



specific mutations. This tenet was challenged by Cairns et al. (1988), who observed that 
under a certain selective pressure, E. coli cells appeared to acquire an excess of beneficial 
mutations. The idea that cells can somehow 'direct' evolution was thought provoking, and 



Bridges 


(1998 


); 


Foster 


(1999 


); 


Cairns 



(1998); Hall (1998); Rosenberg (2001)). While the notion that cells can directly de- 



cide in which genomic regions to increase their mutation rate has been mostly abandoned 



( |Foster[ |1998} |Cairns[|1998[ ), the original observations by [Cairns et aZ.| ( |1988[ ) have been 
corroborated (see above reviews). 

If mutations arise independently of their phenotypic consequences, then how can adapta- 
tions occur that require multiple amino acid mutations and for which the intermediate stages 
are either selectively neutral or disadvantageous? Large populations can climb multiple fit- 



ness peaks, even with disadvantageous intermediate alleles (Behe and Snoke, 2004; Wein 



reich and Chao 2005). Although no new mechanisms are therefore required to explain 



the evolution of complex proteins (Lynch, 2005), we propose that errors in transcription 



and translation (phenotypic mutations) allow the selection of the intermediate mutations of 
a multiple-mutation requiring trait, and can thus speed up the evolution of complex traits. 
Studies on the phenotypic mutation rate indicate that it is orders of magnitude larger 



than the genotypic mutation rate (Springgate and Loeb , 1975 Edelmann and Gal 



lant 



1977): the global phenotypic non-synonymous mutation rate has been estimated to 



be 4.5 x 10~ 4 mistranslations per codon (Ellis and Gallant, 1982; Shaw et al., 2002) 



compared with a genotypic mutation rate of between ~ 10~ 7 to 10 -11 (Drake et al., 1998). 
Consequently, for a protein of 300 residues, on average more than 1 in 10 copies of the protein 
will contain a mutation. Using mutation rates derived from the literature and conservative 
biological assumptions, we show via mathematical modeling and simulations that pheno- 



typic mutations allow evolution to select for neutral intermediate alleles of a multi-mutation 
trait, actually selecting for proteins whose exact DNA sequence is not in the organism under 
selection. Evolution is then able to look ahead for evolutionary jackpots in sequence space. 

Our theory is based on the following hypothetical scenario. A protein can increase the 
fitness of an individual if it evolved a specific trait. This trait requires two mutations, for 
example a disulphide bridge between two cysteine residues. Having only one of the required 
mutations is either selectively neutral or deleterious, however when an individual has only 
one mutation, small amounts of the protein with both mutations will be produced due to 
phenotypic mutations. If the presence of both mutations at low concentrations provides 
even a small fitness improvement then the allele with one mutation will spread though 
the population. As the frequency of the intermediate allele increases, there is a greater 
probability that if the second mutation occurs, it will be the presence of the first mutation, 
and thus provide the full fitness benefit. 

The aim of this article is to show that adaptive phenotypic mutations can undergo positive 
selection under biologically plausible conditions, allowing proteins greater access to features 
involving multiple mutations. We thus introduce the notion of a "look-ahead" effect. The 
name indicates that seemingly unsurmountable evolutionary barriers can be overcome thanks 
to phenotypic mutations which are not yet present in the genome. We wish to emphasize 
that the look-ahead concept is firmly grounded on the idea of chance and necessity and by 
no means insinuates a teleological feature of molecular evolution. 

2 Materials and Methods 

Model assumptions: We model the scenario of a protein evolving a trait that requires 
two mutations. The model is based on a population-genetics framework where a single gene 
can evolve into different alleles. We do not consider duplication and divergence of genes. 
In addition, the process described here will likely only occur for proteins with sufficiently 
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Variables used in this work 



(3 Phenotypic mutation rate from allele 1 to allele 2 

\i Null mutation rate per residue 

N Population size 

r Number of residues available for the second mutation 

s Selection coefficient for allele 2 

U Generational (DNA) mutation rate from allele 1 to allele 2 



long half-lives, as the protein must persist for some time to exert a phenotypic effect. As we 

model only a single gene, we expect our results to be more relevant for single-celled organisms 

and viruses than for multicellular organisms, which tend to have larger genomes and smaller 

effective population sizes than microorganisms. 

Allele Fitness without Fitness with 

(mutations) phenotypic errors phenotypic errors 



_ — 1 1 



1 — 



— 1 1 + sp 



2 — 



— 1 + s 1 + s 



Figure 1: The three alleles (or genotypes). The vertical lines in the genes indicate the 
number of key mutations required for the novel two-residue function. The fitness of the 
allele 1 increases if phenotypic mutations are taken into consideration. 

The model consists of the evolution of three non-recombining haploid genotypes, where 
each genotype contains one of the three alleles shown in Figure [TJ The three different alleles 
are named according to number of relevant mutations, corresponding to zero mutations 
(allele 0), a single mutation (allele 1), and both mutations (allele 2) required for the adaptive 
feature. Having both mutations of the adaptive feature provides a selective advantage s. We 
assume that the intermediate allele (allele 1) is selectively neutral if transcribed and translated 
without error. We specifically take into consideration errors in transcription and translation, 
that is, phenotypic mutations. 

In the model, the population initially consists of one individual carrying allele 1 and 
N — 1 individuals carrying allele 0. So long as allele 1 is present, allele 2 can be generated 
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by mutations. The population evolves for a fixed time period, during which allele 2 can be 
generated by mutation and go to fixation. 

In each generation, selection increases the frequency of the alleles according to their 
corresponding fitness values. Allele has a fitness of 1. Allele 2 has a fitness of 1 + s, 
where s is the selection coefficient provided by the adaptive feature. The fitness of allele 
1, the intermediate allele with only a single mutation, depends on the phenotypic error 
rate. Most phenotypic errors will be neutral or deleterious, however some will be beneficial. 
For simplicity, we assume that the length of the protein and the expression level are both 
constant. In addition, we also assume that the cost due to deleterious phenotypic errors 
is also constant. The effect of these parameters on this model will be the subject of future 
work. If there are no phenotypic mutations, allele 1 has the same fitness as allele 0. However, 
if phenotypic mutations occur, allele 1 can produce a small number of allele 2 proteins due 
to phenotypic errors. The fitness of allele 1 is therefore dependent on the number of such 
errors. 

We assume that fitness is a linear sum of individual proteins, meaning that if some 
phenotypic variants of a protein have a higher fitness, then the overall fitness of that allele 
is proportionally increased. 

We let r be the number of residues that can potentially complement the first mutation 
to provide the full two-residue adaptive feature. These r residues represent, e.g., the sites at 
which the second cysteine of a cysteine bridge could arise; other similar two-residue mutations 
that significantly improve functionality can be proposed. Two residues that comprise an 
adaptive trait are likely to co-evolve, because if a mutation occurs in one of the residues, 
selection strongly favors a compensatory mutation in the other. Based on a large data set, 



Martin et al. (2005) found that co-evolving residues are spatially near. Co-evolving residues 



were, on average, 98.6 amino-acids apart along the sequence, but had a mean spatial distance 
of 6.9 A. This spatial distance can be compared to the width of the van der Waals volume 
of an amino-acid (5-6 A), showing that most co-evolving residues are effectively in contact 
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proximity. Therefore, r is mostly independent of the size of the protein, as long as the protein 



is of sufficient length. Bloom et al. (2006) calculated the mean contact density (the mean 
number of residues in contact with a given residue) for 194 yeast proteins, and found that 
most residues have a mean contact density of seven to eight residues. In this work we use 
r = 8. 

Given r possible positions for the second residue, and assuming that each position requires 
a specific residue, the fraction of proteins of allele 1 containing the second (now highly 
beneficial) mutation is (3 = jgX, where A is the per codon non-synonymous phenotypic 



mutation rate. In this model, we use A = 4.5 x 10 4 mistranslations per codon (Ellis and 



Gallant 


1982 


Shaw et al, 


2002) 



fitness, giving allele 1 a fitness of 1 + s/3. 

When considering genetic (i.e. inherited) mutations, for simplicity we neglect back mu- 
tations (e.g. from allele 1 to allele 0), and assume there are no recurrent mutations of allele 
1 from allele (the model starts with a single copy of allele 1). Allele 2 arises via a mutation 
from allele 1. We ignore the possibility of a double mutation directly from allele to allele 
2, as this probability is extremely small. The genetic mutation rate for allele 1 mutating 
into allele 2 is derived as follows: For microbes, the rate of mutations per nucleotide per 



generation is between ~ 10 7 to 10 11 (Drake et al., 1998). Here we use 10 8 as the 
non-synonymous mutation rate per codon per generation. The resulting mutation rate for 



changing allele 1 into allele 2 is U — ^10 



19 



10" 



Genes can also acquire null mutations, rendering the gene non-functional and therefore 
eliminating the organism. The null mutation rate for protein-encoding genes is on the order 



of 10 per generation (Drake et al, 1998). However, this rate will depend on the length 
(L) of the protein. Assuming an average protein length of 300 residues, the per-residue null 
mutation rate is given by 10 _6 /300 =~ 3.3~ 9 . For a protein of length L, the null mutation 
rate is given by /i = 3.3" 9 L. 
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Simulations: The numerical simulations were written in Java using the Colt scientific 



library (COLT PROJECT, 2007) for the generation of random numbers. The analytic expres- 
sions were evaluated using both Mathematica and Python, the latter in conjunction with 



the SciPy package (SciPy, 2007). Source code for the numerical simulations is available on 
request from DJW. 

The population in each simulation is represented by three numbers, corresponding to the 
abundance of each of the three alleles. As described, the initial abundances are N — 1,1,0 
for alleles 0,1,2, respectively. The simulation runs for a specified number of generations 
T. We used T = 5 x 10 5 throughout this work. Strictly speaking, T is the number of 
generations in which allele 1 can mutate into allele 2; for later generations this possibility 
of mutation is disabled. If allele 2 is present at time T, then the simulation is continued 
until allele 2 is either lost or has reached fixation. Generations are discrete, with mutations, 
selection, and drift occurring at each generation. During each generation we perform the 
following steps. First we check if either allele or allele 2 has reached fixation; if so, we stop 
the simulation, as both cases are absorbing states. Next, for each allele we check for null 
mutations by drawing a random number from the Poisson distribution where the expected 
number of events is the null mutation rate \x multiplied by the total number of individuals 
with the given allele. Mutations from allele 1 to allele 2 are computed in a similar manner, 
where the expected number of events is U multiplied by the number of allele 1 individuals. 
Then, after the possible production of the mutant allele 2, selection acts on the fitness of 
the alleles, where the frequency of each allele is multiplied by its corresponding fitness, 
[1,1 + s/3,1 + s] for alleles [0,1,2], giving the new number of alleles in a possibly larger 
population. Finally, the next population of N individuals is chosen by recursively sampling 
from the binomial distribution, representing random genetic drift. Allele is first sampled 
with the mean= (frequency of allele 0), and the (number of trials) —N. Allele 1 is then 
sampled from the combined allele 1 and 2 individuals. The number of simulations where 
allele 2 becomes fixed is divided by the total number of simulations, giving an estimate of 
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the fixation probability. The number of simulations for each parameter set was between 10 
and 10 9 . 



3 Results 

3.1 Analytical fixation rate of allele 2 

To calculate the fixation rate of allele 2 we have to consider the two fates of allele 1. Firstly, 
allele 1 can become lost. In this case allele 2 can only be generated during the period of drift 
of allele 1. The alternative fate of allele 1 is fixation. Then allele 2 can be generated either 
while allele 1 drifts to fixation or after allele 1 is already fixed. We would like to know how 
many mutation events from allele 1 to allele 2 are expected for either fate of allele 1. We let 
n(s/3) be the expected number of mutation events for when allele 1 is eventually fixed, and 
nioss(s/3) be the expected number of mutation events for the case when allele 1 is lost. We 
can calculate n(s/3) and rii oss (sf3) from diffusion theory, by integrating over the sojourn times 
of allele 1. The corresponding calculations are cumbersome but straightforward, and for the 



sake of brevity we present the details in the Appendix (A. 4 and A. 5). For n(s(3), allele 2 can 
be generated as allele 1 drifts to fixation, and also after allele 1 has already reached fixation. 
For ni oss (s/3), allele 2 can only be generated while allele 1 drifts. 

Assuming that m is the expected number of times allele 2 is generated, what is the 
probability that at least one copy goes to fixation? The probability of fixation of a single 



copy of allele 2 is u(s) KlMURA (1962). (In Appendix A.l , we reproduce the exact expression 
for u(s), as well as approximations for large and small s.) Thus, if allele 2 is generated k 
times, its probability of fixation is 1 — [1 — u(s)] k . Since the probability that allele 2 is 
generated k times follows a Poisson distribution with mean m, we find for the probability v 
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that at least one of the mutations to allele 2 goes to fixation 



k 

1 - e- mu(s) . (1) 



' ' ' kl 

k 



We calculate this probability separately for n(s0) and rii oss (s/3), setting m equal to either 
of these values. We assume that T is sufficiently large so that allele 1 has time to reach 
fixation within this interval (we assume T > 2N). Then the probability U2{s,(3) that allele 
2 is generated and goes to fixation (starting with a single copy of allele 1) is 

u 2 (s,P) = u{s(3)(l - e-^'^W) + (1 - u(s/3)) (l - e ~ ni ^>^). (2) 

The first half of the equation stems from the case when allele 1 eventually reaches fixation, 
where the probability that allele 1 becomes fixed, u(s/3), is multiplied by the probability v 
that at least one copy of allele 2 is generated and fixed. The second half corresponds to 
the case of loss of allele 1 from the population, where the probability of loss of allele 1, 
(1 — u(sj3)), is multiplied by the probability of at least one mutation from allele 1 to allele 
2 and subsequent fixation of allele 2. Taking into account allele 2 mutations during allele 1 
loss is important especially for small s. Allele 1 is more likely to be lost than fixed for small 
s, but can occasionally drift for long times before being lost. 



In the limit f3 — > 0, i.e., in the absence of phenotypic mutations, we find with Eqs. (A2) 



(A27), and (A35) 



U2 ( Sj 0) = ^1 - e -NU(T-N)u(s) /N _ e -NUu(s)_ (3) 

(We assume that N ^> 1, and neglect corrections of order 1 compared to N. Note that we 
cannot simplify (N + l)/N to 1, because for small U, l-e~ NUu ^ and (l- e - NU ^- N ^)/N 
are of the same order in N.) 
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As we are interested in the effect of phenotypic mutations {(3 > 0) compared to the case 
without phenotypic mutations ((3 = 0), we define the increase in the probability of fixation 
from advantagous phenotypic mutations (the look-ahead effect) as 

We can broaden the assumption of T > 2 N to T — > oo with good accuracy. For T — > oo, 
if allele 1 is destined to reach fixation, then the probability of generating at least one copy 
of allele 2 that goes to fixation approaches 1. Therefore, 1 — e~ n ^ u ^ — ► 1, in this limit, 
and thus 

_ u{s(3) + (1 - u(«fl) (1 - e - ni °^ sl3)u{s) ) 
? w (TV + l)/iV - e -w«W ' (5) 

Apart from a correction for the case when allele 2 occurs while allele 1 is destined for 
extinction, Equation (|5| is just the ratio of the probability of allele-1 fixation in the presence 
and absence of phenotypic mutations, u(s/3)/u(0) = Nu(s(3). 
To first order in s/3, Eq. ^ simplifies to (Appendix A. 6) 



^l + Nsf3 + 0(s 2 (3 2 ). (6) 

We can see from this equation that the look-ahead effect becomes important when N is on 
the order of l/(s(3). 

For Ns(3 ^> 1, only the first term contributes to the numerator in Eq. (pi), and we obtain 



(Appendix A. 7) 



(N + 1)/N -exp[-NU(l - e-' 2s )] 



3.2 Simulations 



We confirmed our analytic results for the fixation probabilities U2{s,P) and U2(s, 0) by nu- 
merical simulation, for different values of s (Figure^)). With a population size N = 10 4 , 
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Figure 2: Fixation probability of allele 2 (w 2 ) vs. the selection coefficient s. Black is for 
u 2 (s,P), grey is for u 2 (s,0). Solid lines are predictions according to Eq. (|2| and (|3]), data 
points are for simulations with 10 9 repeats. N = 10 4 , U = ^lO" 8 , (3 = 0.00019, T = 5 x 10 5 . 
Error bars are standard errors. 

the effect of phenotypic mutations can be seen for s > 0.1, and increases for larger s. For 
s < 0.1, the effect is too small and the intermediate allele is effectively neutral, meaning the 
fixation of allele 2 depends on the random fixation of the neutral allele 1. The look-ahead 
effect, £, shows the simulation results compared to Equations (pi), (j6]) and (jT]). Figure [3] 
shows the magnitude of the look-ahead effect for the same parameter settings. For large s, 
the look-ahead effect can inflate the probability of fixation of allele 2 by several orders of 
magnitude. 

We also display the different analytic expressions for £ in Figure [3j The approximation 
(p|, derived in the limit T — > oo, works well for all values of s. The approximation ([6]), derived 
for small s/3, captures correctly the magnitude of s at which the look-ahead effect starts to 
operate, i.e., s > l/(N(3). Similarly, approximation (J7|), valid for Ns(3 ^> 1, approximates £ 
well for larger s. 

Figure [4] shows £ for different population sizes. As expected from the condition s > 
l/(N/3), the look-ahead effect will work with smaller selection coefficients s in larger popu- 
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Selection coefficient s 

Figure 3: Look-ahead effect (£) due to phenotypic mutations vs. the selection coefficient s. 
The solid line is for Eq. (|5]), dashes are for Eq. ([6]), dots are for Eq. QTh , and data points are 
for simulations with 10 9 repeats. N = 10 4 , U = ^10" 8 , (3 = 0.00019, T = 5 x 10 5 . Error 
bars are standard errors. 

lations. For large s, £ saturates at approximately N. 

We studied the effect of different values of the phenotypic error rate (3 (Fig. [5]). As the 
error rate (3 increases, the look-ahead effect £ increases by the same order of magnitude. For 
a very high phenotypic error rate of (3 = 0.019, the look-ahead effect is present for very small 
values of s. However, such a high error rate is likely to be severely detrimental, and in our 
model we do not take into account the loss of overall fitness for increasing phenotypic error 
rates. Conversely for smaller (3, the look-ahead effect is restricted to large s. 
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Figure 4: Look-ahead effect (£) due to phenotypic mutations vs. the selection coefficient 
s for different population sizes (N). Solid lines are from Equation (Jsj) , data points are for 
" repeats. U = 4l0~ 8 , (3 = 0.00019, T = 5 x 10 5 . Error bars are standard 
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Figure 5: Look-ahead effect (£) due to phenotypic mutations vs. the selection coefficient s 
for different phenotypic error rates (/?). Solid lines are from Equation (|5]), data points are for 
simulations with 10 8 repeats. N = 10 4 , U 
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4 Discussion 



We have described a model demonstrating the consequences of positive phenotypic mutations 
on the evolution of a single gene. We have compared numerical simulations with the analyt- 
ical approximations and found them to be in good agreement. When phenotypic mutations 
exert an effect on fitness, selection can operate on the intermediate allele of a complex trait, 
which otherwise (without phenotypic mutations) would be neutral. We refer to selection for 
the intermediate allele as the look-ahead effect, because this effect allows evolution to select 
for sequences not yet in the genome. 

The approximation for small s/3, Eq.([6]), shows most clearly the relationship between the 
parameters. The look-ahead effect is proportional to N, s, and f3, and sets in when N is on 
the order of l/(s(3). For large Ns(3, the look-ahead effect saturates. The asymptotic value 
of £ is approximately N for NU <C 1. Therefore, large populations have two advantages 
over small populations in terms of the look-ahead effect: the effect sets in for smaller values 
of s, and saturates at a larger asymptotic value £. Of course, even in the absence of the 
look-ahead effect, larger populations can more easily traverse multiple local fitness peaks 



( |Weinr.eich and Chaq[ [2005]) . 

Because the selection coefficient s depends on the environment, a valid question is how 
often does s reach sufficiently high levels so that the look-ahead effect can operate. For 
microbial species such as bacteria, sufficiently large s should be reasonably common. Many 



bacteria experience highly fluctuating (Smit et al, 2001) and structured (Baquero and 



Negri, 1997) environments, where growth is limited by the lack of a key trait. An obvi- 



ous and extreme example is antibiotic resistance. Evolving a defense against an antibiotic 



molecule can involve only a few amino acids (Palzkill et al., 1994), and those individuals 



that can generate an enzyme capable of degrading the antibiotic, even if briefly or weakly, 
will have a very large fitness increase. In fact, if the efficacy of the antibiotic is 100% on 
susceptible genotypes, a mutation providing only moderate resistance has an infinite selec- 
tive advantage. And even for very small antibiotic concentrations, mutants differing by only 
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two amino acids at a single /3-lactamase gene can be selected effectively (Baquero et al. 



1997, 1998). Thus, bacteria may frequently experience environments in which a large fit- 



ness increase (large s) is only a few mutations away. Similarly, in bacteriophages, selective 



coefficients s of 10 or more are not uncommon, even for individual mutations (Bull et al. 



2000). 



Our work is entirely theoretical, but we expect that it will be possible to experimentally 
verify our predictions in future work. For experimentally observing the look-ahead effect, 
we would need a system where s and N are both large, while /3 (the phenotypic mutation 
rate) can be modified. The values of both N and s used in this work are well within bio- 
logically realistic ranges achievable in a microbiological laboratory. Conditions for large s 
may be created with e.g. antibiotic resistance, which is a common laboratory workhorse. 



Unfortunately, many antibiotics function by reducing translation fidelity (Ogle and Ra- 



MAKRISHNAN 



2005), and thus would conflate s and (3. Changing (3 could involve a mutated 



ribosome. Ribosomes appear to be optimized for accurate and efficient translation of mRNA 



( Baxter- Roshek et al. , 2007), and several examples of altered ribosome fidelity exist, both 



increasing ( |Vila-Sanjurjo et al] |2003[ ) and decreasing fidelity ( |Friedman et al] [1968J). 
Specific regions of the ribosome rRNA sequence have been identified as influencing fidelity 



(O'Connor et al, 1997), and various agents can reduce fidelity, e.g., streptomycin, mag- 



nesium, and ethanol (Friedman et al, 1968). Few mutations may be sufficient to alter 
the fidelity of a ribosome, for example, a single mutation in the S5 ribosomal protein in E. 



coli increases frameshifting and nonsense mutations (Kirthi et al., 2006). In yeast, muta- 
tions in the 18S RNA have been found that both increase and decrease translational fidelity 



( |Konstantinidis et al] [2006J). 

In this work, we have calculated the look-ahead effect from a comparison between the 
two cases of (3 > and (3 = 0. The latter may not be experimentally possible; any exper- 
iment will likely compare two different positive values of f3. Nevertheless, Figure [5] shows 
that a larger look-ahead effect can be achieved with a higher /3, where increasing (3 by one 
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order of magnitude both increases the look-ahead effect by an order of magnitude and low- 
ers the smallest s where an effect is observed. Of course, our model does not take into 
account the loss of fitness or other confounding effects from a higher phenotypic mutation 
rate. Thus, a balance must be found in having two different values of (3 that are different 
enough to measure, while at the same time minimizing the confounding effects. The most 
obvious consequence of increasing the phenotypic mutation rate is that overall fitness may 
be reduced, for example in E. coli, where a higher translational error rate activates stress 



responses (Fredriksson et al. , 2007), or in mouse, where such errors are implicated in neu- 



rodegeneration caused by misfolded proteins that aggregate (Lee et al., 2006). Increasing 
translational fidelity may not come without fitness cost either. The hyperaccurate mutations 



in the 18S RNA in yeast (Konstantinidis et al. , 2006) cause an increase in oxidative stress. 
This observation suggests that cells consume more energy to achieve hyper accuracy. It may 
also partially explain why the phenotypic error rate is much higher than the genotypic error 
rate, as there is possibly a direct disadvantage in reducing the phenotypic error rate, rather 
than only reducing the selective advantage that occurs if the phenotypic error rate is reduced, 



as discussed in Buerger et al. (2006). 



Buerger et al. (2006) asked whether evolution has selected for the current phenotypic 



error rate, which does not differ significantly between eukaryotes and prokaryotes (ELLIS 



and Gallant, 1982 Shaw et al., 2002) even though the source of errors is different. They 



suggested that the increase in fitness becomes incrementally smaller for improvements to 
transcription and translation fidelity. We would like to speculate that the phenotypic error 
rate is on the border between minimal costs (of e.g. misfolded proteins) and maximum payoff 
(via the look-ahead effect). 

The goal of our analysis was to demonstrate that the look-ahead effect is theoretically 
possible, and as such, we intentionally excluded confounding factors for the sake of clarity. 
There are several aspects not considered in our model that may play important roles. For 
example, in this work we did not consider the expression level. For low expressed genes, the 
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mutation from allele 1 to allele 2 will occur less frequently compared to highly expressed 
genes. However, if allele 2 is produced it will be at a higher concentration (of allele 2 mutant 
proteins in a population of allele 1 proteins), as the overall copy number of allele 1 is low. 
This difference in expression levels is likely reduced in a large population, where beneficial 
mutations occur with sufficient frequency. Another factor related to the expression level is 
translational robustness. It has been proposed that highly expressed genes are under selection 



to properly fold despite phenotypic mutations, and consequently evolve slower (Drummond 



et al, 2005; Wilke and Drummond, 2006). If a gene is robust to translational errors, 



then it can tolerate a larger variety of mutations, of which some may be intermediates to a 
new adaptive multi-residue trait. Thus, translational robustness may increase the sequences 
available for experimentation at the phenotypic level. However, if the intermediate allele is 
itself not robust to errors in translation, then it will not be neutral, and may be selected 
against. The location of the protein trait will also influence the viability of the intermediate 
allele: mutations near the surface of the protein are less likely to disrupt the protein compared 



to mutations in the core flToKURlKl et al.\ |2007j). 

In conclusion, we propose that organisms can experiment with protein sequences that 
are mutationally close to the current sequence, but not yet in the genome. This effect allows 
selection for intermediates of complex traits, opening up a more direct route to the trait and 
thus reducing the time needed for fixation in the population. 
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A Appendix 

Here, we present the details of our analytic derivations. 



A.l Probability of fixation 



According to (KlMURA, 1962), the probability of fixation u(s) of a single allele with selection 



coefficient s is given by 

1 - e~ 2s 
1 -e 



u ( s ) = i ^=2Nl- ( A1 ) 



For s < 1/N, this expression simplifies to 

u(s) = ^ + ^s + 0(s 2 ), (A2) 
whereas for Ns ^> 1, this expression simplifies to 

u(s) « 1 - e- 2s . (A3) 

A. 2 A single allele drifting to fixation or loss 

We first consider a single allele with selective advantage s drifting to fixation or extinction, 
and ask how many mutations this allele generates until it is either fixed or lost. We will treat 
these two cases separately. Let n fix (s) be the expected number of mutations generated while 
the allele drifts to fixation, and let n\ oss (s) be the expected number of mutations generated 
while the allele drifts to extinction. We calculate these two quantities using diffusion theory, 
by integrating the sojourn times of the allele over all frequencies. 



For an allele with selective coefficient s and starting at frequency p = 1/N, |Nagylaki 



(1974) calculated its mean sojourn time r(y) between frequencies y and y + dy as 



r(y) = 2{V(y)G(y)}- l {u lo M/N)g(0,y)6(l/N - y) + u Rx (l/N)g(y, 1)% - 1/N)]. (A4) 
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Here, 



V(y)G( V )=y(l-y)e- 2N ^/N, (A5) 

e -2Nsa _ e ~2Nsb 

</M) = — , (A6) 

e -2Nsp _ q-^Ns 

Mloss(p) = ; ZZ^l , ( A7 ) 



l-e~ 

l _ £ -2Nsp 

1 -e- 



Wfix(p) = 1 - Wi oss O) = 1 IZ^7> ( A ^ 



and ^(z) is the Heaviside step function. We want to integrate expressions involving r(y) 
from y = to y = 1. Since y = 1/N corresponds to a single copy of the allele that drifts to 
fixation, values of y less than 1/N are not relevant for our analysis. Therefore, we discard 



the term proportional to 6 (1/N — y) in Eq. (A4), and use in what follows 



r(y) = 2ufr(l/N)g(y, l)/[V(y)G(y)] for y > 1/N. (A9) 
A. 3 Number of mutations conditional on fixation 



For the sojourn time conditional on fixation, r^ x (y), Nagylaki (1974) finds 

rsx(y) = r{y)ufa{y)/ufa{p). (A10) 

Using this expression, we have 



n fix (s) = nu I T &x (y)ydy. (Ail) 

Jl/N 



'1/N 

Plugging the expressions for V(y)G(y), g(a,b), Ufix(p), and r(y) into rg x (|/), we arrive at 

1 (1 - e- 2Ns y)(l - e- 2Ns ^) ,. 1n . 
= 7{ T^- ) - ylhj) < A12 > 
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This expression corresponds to the one by Ewens (1973). Note that yr^ x (y) — > for y — > 0. 



Therefore, we can extend the lower limit of integration to in Eq. (All ), and rewrite ng x (s) 
as 



NU 



n Rx (s 



s(l-e- 



2Ns' 



-I(2Ns) 



(A13) 



with 



,. (1 - e -. )(1 - e -^-,») 

Jo 1-2/ 



(A14) 



The integral 1(a) can be rewritten as 



1(a) = 7 - Ei(-a) + ln(a) + e~ a [7 - Ei(a) + ln(a)], 



(A15) 



where 7 ~ 0.5772 is the Euler-Mascheroni constant and Ei(z) is the exponential integral, 



Ei(z) 



00 e - t 



dt. 



(A16) 



For s < 1/iV", we find 



n &x (s) = N 2 U + 0(s 2 ). 



(A17) 



For Ns ^> 1, we obtain the asymptotic expansion 



NU 

n^s) « pn(2i\Ts) + 7], 



(A18) 



using Abramowitz and Stegun (1964) 5.1.51 



Ei(-z) ~ 



1 2 6 

z z z z a 



for large z. 



(A19) 
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A. 4 Number of mutations conditional on extinction 



For the sojourn time conditional on extinction, Ti oss (y), [Nagylaki (1974) finds 

TiossO) = r(y)u loss (y)/u loss (p). (A20) 

Using this expression, we have 

^ioss(s) = NU / n oss (y)ydy. (A21) 

Jl/N 

Plugging the expressions for V(y)G(y), g{a,b), Ui oss (p), and r(y) into r loga {y), we find 

1 e 2s -l ( e ~ 2Nsy - e~ 2Ns ){\ - e - 2Ns 0--v)) 

T ^s(y) = s(1 _ e _ 2JVs) i _ e -2(N-i )s y(l-y) ' (A22) 



We rewrite ni oss as 



with 



NU e 2s - 1 

nioss = s ( 1 _ e -2Ns) 1 _ e -2(N-l)s J ( N > S ) ( A23 ) 



fl ( p -2Nsy _ -2Ns\(-t _ -2Ns(l-y)\ 

W s) = / * ' K Uy. (A24) 

Ji/n l ~y 

The integral can be rewritten as 

J(N, s) = -2e- 2Ns (7 - Chi[2(A - l)s] + ln[2(A - l)s]) , (A25) 
where Chi(z) is the hyperbolic cosine integral, 

Chi(z) = 7 + + T C ° sh(t) ~ 1 dt. (A26) 



For s < l/N, we find 

n loss ( S ) = (A-l)f/ + 0( S 2 ). (A27) 
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For Ns ^> 1, we obtain the asymptotic expansion 

U 



2s 



using 



™ioss(s)~— (l-e^ 2s ), (A28) 



Ei(z) e z 

Chi(z) « — « — for large 2. (A29) 
[This expansion follows directly from the definitions of Chi(z), cosh(z), and Ei(z).] 



A. 5 Number of mutations within a given time interval 

We now extend the derivations in Section IA.3I to calculate the number of mutations to allele 
2 generated within a certain time interval T, conditional on fixation of allele 1. We assume 
that T is sufficiently large so that allele 1 has time to reach fixation within this interval. We 
only consider the case conditional on fixation because no new mutations are generated once 
allele 1 has gone extinct. 

We calculate n(s) = na x (s) + nx(s), where n-r(s) is the total number of mutations gen- 
erated once the first mutation has reached fixation. We have 

n T (s) = NU[T-t fix (s)}, (A30) 

where ts x (s) is the time to fixation of a mutation with selective advantage s. This time is 
given by the integral over all sojourn times, 



1 h(2N s) 
sfl-e- 



Ms) = / Tz x (y)dy = ,/ v _J (A31) 



with 



h{a) = I U j 1 g Uy. (A32) 







A partial fraction decomposition of the integrand reveals that /2(a) = 21(a), and thus we 
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have 



2I{2Ns) 



s(l -e 



Combining this result with Eqs. (A13) and (A30), we find 



n(s) = n fix (s) + n T (s) = NU 
= NUT-n &x (s). 



T 



I(2Ns) 



s(l - e- 2Ns ) 



Note that n(s) = n &x (s) for T = tfix(s). 
For s < 1/N, we find 

n(s) = NU(T-N) + 0(s 2 ). 



For As ^> 1, using Eqs. (A15) and (A19), we obtain the asymptotic expansion 



ln(2As) + 7 



n(s) w NU(T- 

\ s 



A.6 £ for s/? < 1 



From Eq. (El), using Eqs. (A27), (A35), and (A2), we obtain to first order in s(3 



e -NUu(s) _ e -NU(T-N)u(s) 



2n2\ 



u 2 (s,0) 



If further NU(T - N)u(s) < 1, we obtain 



£ « 1 + N(l - jV/2>/3 + 0(s 2 p 2 ), 



and for T — ► oo, 



l + A^s/3 + C(s^). 
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A.7 f for Ns(3 > 1 



For Ns(3 3> 1, only the first term contributes to Eq. (|2]), and we obtain from Eqs. (A36) and 



(A3) 



«aM) = (l 



-2s P\ 



1 - exp ( - NU 



T 



ln{2Ns(3) + 7 

s(3 



(1 - e"-) 



(A40) 



Likewise, in this limit we can simplify Eq. ^ to 



u 2 (s,0) 



N + l 

N 



exp[-NU(T - N)(l - e~ 2s )}/N - exp[-NU(l - e~ 2a )\, (A41) 



giving 



'l-e 



-23(3} 



1 - exp - NU 



T 



\n(2Ns[3) + 7 

s/3 



(1 - e-») 



(N + 1)/N - exp[-NU(T - N)(l - e~ Zs )}/N - exp[-NU(l - e" is )] 



(A42) 



Furthermore, for T — > oo, this expression simplifies to 



(1 - e~ 2s ?) 



(N + 1)/N - exp[-NU(l - e~ 2s )] 



(A43) 



If NU < 1, then £ — > AMn the limit s -> oo. 
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